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Abstract. 

Dynamical decoupling (DD) is a technique for preserving the coherence of quantum 
mechanical states in the presence of a noisy environment. It uses sequences of 
inversion pulses to suppress the environmental perturbations by periodically refocusing 
them. It has been shown that different sequences of inversion pulses show vastly 
different performance, in particular also concerning the correction of experimental pulse 
imperfections. Here, we investigate specifically the role of time-reversal symmetry in 
the building-blocks of the pulse sequence. We show that using time symmetric building 
blocks often improves the performance of the sequence compared to sequences formed 
by time asymmetric building blocks. Using quantum state tomography of the echoes 
generated by the sequences, we analyze the mechanisms that lead to loss of fidelity 
and show how they can be compensated by suitable concatenation of symmetry-related 
blocks of decoupling pulses. 
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1. Introduction 

Dynamical decoupling (DD) [H |2] is becoming an established technique for preserving 
quantum states from decoherence with possible applications in quantum information 
[a H El El [7] and magnetic resonance [H El [TOl El [121 [131 [H] . The technique was devised 
to increase decoherence times by refocusing the system-environment interactions with 
a sequence of control pulses periodically applied to the system. Recent experiments 
have successfully implemented DD methods and demonstrated the resulting increase of 
coherence times in different systems [l5l[l6l[I71[l8l[T9l[20l[2ll[22l[23]. These works also 
showed that the performance of DD sequences can be limited or even counterproductive 
if the accumulated effect of pulse imperfections becomes too strong [181 [20l [22l [23] . Oiig 
approach to compensate the effect of these errors is to combine one basic decoupling 
cycle with a symmetry- related copy into a longer cycle. The resulting cycle can be more 
robust, i.e. less susceptible to pulse imperfections than its building blocks, provided the 
basic blocks are well chosen and combined in a suitable way. 

In the field of nuclear magnetic resonance (NMR), symmetry- related arguments 
have often been used for constructing supercycles [241 [2Sl [26l [27]. Using the symmetry 
properties of specific interactions, it is possible to remove them selectively while retaining 
or restoring others [281 [29] . Symmetrization is widely used to eliminate unwanted odd- 
order average Hamiltonian terms [30]. This approach has been instrumental in the design 
of high-performance decoupling and recoupling sequences [281 [22]. Besides sequence 
development, symmetry arguments have also been used extensively in the design of 
individual pulses with reduced sensitivity to experimental imperfections [3T1 [32]. 

The main goal of this paper is to investigate differences between otherwise identical 
DD cycles, in which the timing of the pulses is either symmetric with respect to time 
reversal, or not. As the basic block we consider the XY-4 sequence [33]. We compare 
the performance of the basic sequences as well as compensated higher-order sequences 
and analyze their imperfections theoretically by average Hamiltonian theory [231 [M] and 
experimentally by applying quantum state tomography [35l [36] to the system after the 
end of each decoupling cycle. 

This paper is organized as follows. In Section [2] we introduce the basic idea of 
dynamic decoupling and demonstrate the relevance of time reversal symmetry in this 
context. In the subsequent sections [31 and [H we compare different sequences based on 
symmetric or asymmetric building blocks. In the last section we draw some conclusions. 

2. Symmetrization in DD 

Dynamical decoupling is a technique in which the coherence of qubits is dynamically 
protected by refocusing the qubit-environment interaction [Il[2]. Within this technique, 
a sequence of vr rotations is periodically applied to the system. For a purely dephasing 
environment, i.e. one that couples only to the 2;-component of the system qubit, this can 
be achieved simply by a train of identical vr-pulses, the so-called Carr-Purcell (CP)- [37] 
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or Carr-Purcell-Meiboom-Gill (CPMG)-sequence [38] • The shortest DD sequence that 
cancels the zero-order average Hamiltonian for a general system- environment interaction 
[11 [39] is the XY-4 sequence [33] (see figures [TH and Wp)- This sequence also has the 
advantage of being much less sensitive to pulse imperfections than the CP-sequence 

[HHlilo]. 
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Figure 1. Schematic representation of dynamical decoupling sequences: (a) time 
symmetric XY-4 and (b) asymmetric XY-4. 

In the spectroscopy and quantum computing communities, two versions of the XY- 
4 sequence are used that differ by a seemingly minor detail. As shown in Fig. [1^, the 
basic cycle originally introduced in NMR [33l|40] starts with a delay of duration r/2 and 
ends with another delay of the same duration. It therefore shows reflection symmetry 
in the time domain with respect to the center of the cycle. In contrast to that, the 
sequence used in the quantum information community [H [391 El il2l US] starts with a 
delay of duration r and ends after the fourth pulse (see Fig. [lb). Clearly, this cycle 
is not symmetric in time. One consequence of this small difference is that in the case 
shown in Fig. [1^, the echoes are formed in the center of the windows between any two 
pulses, while in Fig. [T]d, the echoes coincide with every second pulse. The separation in 
time between the echoes is therefore twice as long in this case. 

Figure [2] illustrates this difference with an experimental example. Here, we 
measured the time evolution of the ^^C nuclear spin polarization during a CPMG 
sequence, using in one case a time-symmetric and in the other case a non time-symmetric 
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Figure 2. Experimentally observed evolution of the ^'^C nuclear spin magnetization of 
adamantane during the symmetric (red circles) and asymmetric (black squares) CPMG 
pulse sequences. 



cycle. The sample used for this experiment was polycrystalline adamantane. The 
dephasing of the nuclear spin signal originates from the interaction with an environment 
consisting of nuclear spins. To make this environment appear static and generate a 
long train of echoes, we applied a homonuclear decoupling sequence to the protons [21] . 
As expected, in the symmetric case, the echoes appear with half the separation of the 
asymmetric case. 
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Figure 3. ^'^C nuclear spin magnetization after one cycle of symmetric and asymmetric 
CPMG sequences for different cycle times 2 r. 



The larger separation of the echoes also can lead to a faster decay of the echo 
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Table 1. Dynamical decoupling sequences. The top line shows the time symmetric 
and asymmetric versions of XY-4, which can be used as building blocks for other 
sequences. X and Y represent 7r-pulses around the x and y axes respectively. [/"^ is 
the time-reversed sequence and U stands for the sequence with inverted phases. 

amplitude if the environment is not static. This is illustrated in Fig. El which shows 
the decay of the echo amplitude as a function of time. In this case, we did not apply 
homonuclear decoupling and the system-environment interaction is therefore modulated 
by the homonuclear magnetic dipole-dipole interaction between the environmental 
protons |2Tj. The plotted signal shows the echo amplitude measured after a single 
CPMG cycle for the symmetric and asymmetric function of the cycle time 2r. 

Clearly the symmetric cycle is more efficient in preserving the state of the system - in 
agreement with findings from multiple pulse NMR [28l [29] . 

Since any multiple pulse cycle suffers from imperfections and non-ideal properties, 
it is often desirable to construct longer cycles that have better properties than simply 
repeating the basic cycle. Examples of DD sequences that can be constructed from the 
XY-4 cycle include the XY-8 and XY-16 |10] sequences shown in Tabled] Here, the XY- 
8 sequence concatenates an XY-4 cycle with its time- reversed version [40l 1411 l42l |43] , 
thus generating a new cycle, which is inherently time- symmetric, independent of which 
version of the XY-4 sequence was used for the building blocks. 

In the following sections [3] and |4] we show that two symmetric sequences, constructed 
according to the same rules from a basic XY-4 block can have different behaviors 
depending if the basic block is chosen to be symmetric or not. 

3. Average Hamiltonian Theory 

In this section we compare DD sequences constructed from time symmetric and 
asymmetric building blocks in the framework of average Hamiltonian theory [211 [M] . 
While the zero-order average Hamiltonian of the asymmetric sequence is the same as 
that of the symmetric sequence, this is no longer the case for the higher order terms. In 
particular, if a cycle is symmetric, that is if H{t) = Hijc — t) for < t < Tc, where H{t) 
represents the Hamiltonian in the toggling frame [241 El], it can be shown that all odd 
order terms of the average Hamiltonian vanish [30] . Clearly, this condition can only be 
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fulfilled, if the timing of the sequence is symmetric, as in the example of Fig. [T^. 

We first consider the sequence XY-4, which is our basic building block. Our system 
consists of a single qubit, which we describe as a spin 1/2, and an environment, which 
consists of a spin-bath. The Hamiltonian describing the system plus environment is then 

H = Hs + HsE + He, (1) 

where Hs = ujsSz is the system Hamiltonian, S = {Sx, Sy, Sz) is the spin vector operator 
of the system qubit and us is the Zeeman frequency of the system. He is the environment 
Hamiltonian, which does not commute with Hse in general but is not specified further. 
HsE is the system-environment interaction: 

HsE=Y.bkSzl'z- (2) 

k 

In the following, we work in a resonant rotating frame, where Us = and therefore 
Hs = 0. I**^ = {ll^,Iy,I^) is the spin vector operator of the fc*^ environment spin, 6^ is 
the coupling constant between the system and the /c*^ spin of the environment. 

In our case, the dominant source of experimental imperfections are flip-angle errors. 
The actual pulse propagator for a nominal vr rotation around an axis defined by is 
then 

R((j)) = e-^(i+^)-^^ (3) 

where e is the relative flip angle error, 5*0 = cos(f>Sx + sincpSy, and is the phase of 
the pulse. We can write the zeroth {Hq) and first {Hi) order terms of the average 
Hamiltonian for the time symmetric XY-4 sequence 

h! = He (4) 

H!=^-^Sz-j:b,'^{Sx + Sy)I^. (5) 

Details of the calculation are given in the appendix. The zeroth order average 
Hamiltonian matches exactly the target Hamiltonian, and for perfect pulses (e = 0), 
the first order term vanishes, Hf = 0, as expected for any symmetric sequence. For 
finite pulse errors, the first-order term contains a rotation of the spin qubit around the z 
axis by an angle 5e^7r^/4. This term results from the accumulated flip angle errors and 
is independent of the environment. Since this term is proportional to the square of the 
flip angle error e, it generates a rotation in the same direction for all spins, independent 
of the actual field that they experience. 

The second term in eq. ([5]), in contrast, is linear in e. For an optimal setting of 
the pulse, e is distributed symmetrically around zero and the resulting evolution due 
to this term does not lead to an overall precession, but to a loss of amplitude. This 
term combines pulse errors and the system-bath interaction. It arises from the fact 
that pulses that do not implement a vr rotation cannot properly refocus the system- 
environment interaction. 
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Now we compare these results with the average Hamiltonian of the time asymmetric 
form of XY-4: 

11^ = He (6) 

1 IQ^ ' ^ z 

+ tTS,J2Mlz,HE], (7) 
k 

The most striking difference from the symmetric case is the appearance of a new term 
which is a commutator between the internal Hamiltonian of the environment He and 
the system-environment interaction Hamiltonian Hse- 

Under ideal conditions, the first-order average Hamiltonian vanishes for the 
symmetric building block, but not for the asymmetric case. For the asymmetric case 
the third term, which is proportional to [Hse, He] remains. The commutator describes 
the time dependence of the system-environment interaction due to the environmental 
Hamiltonian He- This difference from the symmetric case may be understood in terms 
of the different positions of the echoes shown in Fig. [T] and Fig. [2l In the asymmetric 
sequence, the time between echoes is twice as long as in the symmetric sequence, 
which means that a time-dependent environment has a bigger effect. In the symmetric 
sequence, the effect of the time-dependent environment appears only in the next-higher 
order term. 

Rules for improving the performance of multiple pulse sequences were discussed, 
e.g., in the context of broadband heteronuclear decoupling [H] or for the compensated 
Carr-Purcell sequences |40]. If we combine a XY-4 cycle with its time- reversed image 
to an XY-8 cycle, we obtain a time-symmetric cycle, even if we start from the non- 
time symmetric block. Nevertheless, we expect different results for the two cases. 
An explicit calculation of the average Hamiltonian for the combined cycle shows that 
Hq = Hq = He and Hf = H^ = 0, i.e. all deviations from the ideal Hamiltonian 
vanish to first order. This remains true for finite pulse errors: the symmetry of the 
sequence cancels the effect of pulse errors in all odd-order average Hamiltonian terms. 

We therefore proceed to calculate the second order terms. For simplicity, we do not 
calculate the general expression, but consider two limiting cases. First, we assume that 
the environmental Hamiltonian vanishes. He = 0. The second order term then becomes 

k 

HA = HS^y2^3lSy. (9) 

Again, the average Hamiltonian for the sequence built from asymmetric blocks contains 
an additional error term, which depends on the pulse error and the square of the system- 
environment interaction. 

As the second limiting case, we assume ideal pulses but non- vanishing environmental 
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Hamiltonian, He 7^ 0. The second order terms then become 



Hi = ^[[He,Hse],He~1-Hse] 
o o 

H^ = Hf+^[[HE, Hse]JHe - Hse]- 



(10) 



(11) 



As for the XY-4 sequence, the time-dependence of the environment, represented 
by the commutator [He, Hse] has the bigger effect if the sequence uses an asymmetric 
building block and therefore generates echoes with bigger time delays between them. 

4. experimental Results 
4-1- Setup and system 

For the experimental tests we used natural abundance ^^C nuclear spins in the CH2 
groups of a polycrystalline adamantane sample as the system qubit. The carbon spins 
are coupled to nearby nuclear spins by heteronuclear magnetic dipole interaction 
corresponding to Hse- The protons are coupled to each other by the homonuclear 
dipolar interaction, which corresponds to He and does not commute with Hse- The 
system environment interaction is therefore not static and the carbon spins experience a 
fluctuating environment [18]. Under our conditions, the interaction between the carbon 
nuclei can be neglected and the decoherence mechanism is a pure dephasing process [T8] . 
the evolution of the system and environment is thus described by the Hamiltonian ([1]). 

The experiments were performed on a homebuilt 300 MHz solid-state NMR 
spectrometer. The basic experimental scheme consisted of: a state preparation period, 
during which we prepared the carbon spins in a superposition state oriented along the 
y direction, a variable evolution period, where DD sequences were applied, and the 
final read out period where we determine the final state by quantum state tomography 
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Figure 4. Decay of the My-magnctization for the symmetric version of XY-4 and 
different delays (r) between pulses. 
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Figure H] shows the signal decay for the symmetric version of XY-4 and different 
delays r between the refocusing pulses. For the shortest cycle times, we observe shorter 
decays and oscillations. As discussed in Ref. |18], this is an indication that in this 
regime pulse imperfections play the dominant role. From the decay curves, we extract 
decay times as the times where the magnetization has decayed to 1/e of the initial value. 



4-2. Measured decay times 

Figure [5] shows the decay times of the M^^-magnetization as a function of the delay r 
between the pulses. For delays between 200 fis and 50 /is, the decoupling performance 
improves as the delays between the pulses are reduced. However, as the delay between 
the pulses becomes shorter than 50 /is, the decay time decreases again, in agreement 
with what was observed in [TSl [22] : in this region, pulse errors become more important 
than the coupling to the environment. This occurs equally for both, the symmetric and 
the asymmetric XY-4 sequence. 

If we concatenate the XY-4 sequence with its time-reversed version to the XY-8 
sequence, we obtain qualitatively different behavior for the two different versions of XY- 
4: If we start from the symmetric form of XY-4, the resulting XY-8(S) sequence shows 
improved decoupling performance for increasing pulse rate, without saturating. This 
is a clear indication that in this case, the concatenation eliminates the effect of pulse 
imperfections and generates a robust, well-compensated sequence. In strong contrast 
to this, concatenation of the asymmetric version of XY-4 to XY-8 (A) does not lead 
to a significant improvement: the decay times for XY-8 (A) are identical to those of 
the two XY-4 sequences, within experimental uncertainty. A further concatenation 
to XY-16 does not change this behavior. The qualitatively different behavior of the 
sequences using symmetric versus asymmetric building blocks clearly shows that for the 
asymmetric versions, the pulse errors dominate, while the symmetric ones compensate 
for the effect of pulse errors. 
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Figure 5. Decay time of the Afj,-magnetization for different XY-n sequences. 
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Figure 6. Evolution of the magnetization during the symmetric (A) and asymmetric 
(S) versions of the XY-4 sequence for pulse spacings of t = lOjUs. The top panel shows 
the evolution of the magnetization components as a function of time. The bottom 
panel represents the Bloch vector in the xy plane at different times. The color code 
in the lower panel denotes the time evolution, blue for the initial state and red for the 
final state. 



For a more detailed picture of the process that reduces the signal for high pulse rates, 
we applied state tomography of the evolving qubit by measuring all three components 
along the x, y and ^-direction. Figure |6] shows the observed data for both versions of the 
XY-4 sequence. The oscillation of the x- and y components and the constantly small 
value of the z-component are a clear indication of a precession around the z-axis, in 
addition to the loss of signal amplitude. This combination of precession and reduction 
of amplitude is also shown in the lower part of Fig. [6], where the arrows show the xy- 
components of the magnetization for different times during the sequence. According to 
eqs (IS]) and (l7l), the precession around the z-axis originates from the pulse error term 
5e^7r^/(16r)S'2, which is proportional to and is the same for the symmetric and the 
asymmetric sequence, in excellent agreement with the observed behavior. 

Figure [7] shows the corresponding data for the two XY-16 sequences. Here, as well 
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Figure 7. Evolution of the magnetization during the symmetric (A) and asymmetric 
(S) versions of the XY-16 sequence for pulse spacings of t = 10/is. The top panel shows 
the evolution of the magnetization components as a function of time. The bottom panel 
represents the Bloch vector in the xy plane at different times. The color code is the 
same as in Fig. [6l 



as in the case of XY-8 (data not shown), we also observe a precession for the XY-16 (A) 
sequence, but for the sequence with symmetric building blocks, the oscillation is not 
observed. Again, these results indicate that the sequences built from symmetric XY-4 
blocks have smaller average Hamiltonians and therefore show better performance than 
those built from asymmetric blocks. 

If we change the spacing between the pulses, the behavior remains the same. In 
Fig. [HI we show the measured precession angle around the ^-axis divided by the number 
of pulses. The precession is indistinguishable from zero for the compensated XY-8(S) 
and XY-16(S) sequences. For other sequences, it is significant and independent of the 
delay between the pulses. 

The precession of the magnetization around the 2;-axis that we observe for some of 
the sequences causes a deviation of the system from the desired evolution and reduces 
therefore the fidelity of the process. However, compared to a dephasing process, it is 
easier to correct and can in principle be compensated if it is known. We therefore 
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compared not only the reduction of the magnetization amphtude along the initial 
direction, but also the total magnetization left in the system, which eliminates the 
effect of the precession. Figure [9] shows the decay times of the total magnetization for 
different XY-n sequences. For short delays between the pulses, the difference between 
sequences built by symmetric and asymmetric building blocks is small, indicating that 
the main difference is related to the precession originating from the pulse errors, which 
is better compensated by concatenating symmetric building blocks. For pulse delays 
longer than r ^ 15/iS, we start to see again that the symmetric versions of XY-8(S) and 
XY-16(S) are superior to the asymmetric versions. At this point, the time dependence of 
the environment plays a bigger role and reduces the efficiency of the refocusing [T8l [22]. 
In agreement with eqs IQ and ffTTl) . these effects are bigger for those sequences that use 
asymmetric building blocks. 
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Figure 8. Average precession angle of the transverse magnetization during the 
different XY-n sequences. 
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Figure 9. Decay time of the total magnetization for XY-n sequences. 
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5. Discussion and Conclusions 

Dynamical decoupling is becoming a standard technique for extending the lifetime 
of quantum mechanical coherence. Many different sequences have been put forward 
for reducing the effect of the environmental noise on the system. Since the number 
of possible sequences is infinite, a relatively straightforward approach for designing 
improved sequences consists in concatenating different building blocks in such a way that 
the resulting cycle has a smaller overall average Hamiltonian than that of its component 
blocks. In this paper we consider the XY-4 sequence as a basic building block. Since 
different versions of the XY-4 sequence were proposed in the literature, some of them 
symmetric in time, others asymmetric, we compare these two versions and in particular 
the different sequences that result when they are concatenated with time-inverted and 
phase-shifted copies. Since time-symmetric sequences generate average Hamiltonians 
in which all odd-order terms vanish for ideal pulses, it is expected that they perform 
better than non-symmetric but otherwise identical sequences. Experimentally, we could 
not verify this for the XY-4 sequence, since pulse errors dominate the behavior under 
our experimental conditions. However, in the case of the CPMG sequence, where pulse 
errors are insignificant, we could clearly verify this expectation. 

The symmetry of the basic building blocks is also important when they are 
concatenated to higher order sequences, such as the XY-8 and XY-16 sequences |40]. In 
this case, the odd order terms vanish in the average Hamiltonians of both sequences, 
but the second order terms of the sequences that are built from asymmetric blocks 
contain additional unwanted terms. The experimental data are in agreement with 
this observation: sequences consisting of time symmetric building blocks perform 
significantly better than the corresponding sequences formed by time asymmetric blocks. 

In order to understand the decay processes during the DD sequences, we performed 
quantum state tomography as a function of time. The results from these measurements 
show two different contributions to the overall fidelity loss: A precession around the z- 
axis, which we could attribute to the combined effect of flip-angle errors and an overall 
reduction of the amplitude, which results from the system-environment interaction. For 
short delays between the pulses and correspondingly large number of pulses, the pulse 
error term is the dominating effect. Again, the symmetric and asymmetric version of 
the XY-4 sequence show similar performance. However, as we use them as building 
blocks of the higher-order XY-8 and XY-16 sequences, we find that the effect of the 
pulse errors is almost perfectly compensated if we use the symmetric building blocks, 
while a significant effect remains when asymmetric blocks are concatenated. 

While we have analyzed the effect of symmetry mostly for the XY-tt, sequences, this 
can clearly be generalized. As we showed in Fig. [3l the symmetric version of the CPMG 
sequence shows significantly better decoupling performance than the asymmetric version. 
The same concept can also be applied to the CDD sequences, which are generated by 
inserting XY-4 sequences inside the delays of a lower-order CDD sequence [39]. The 
conventional concatenation scheme [39] uses asymmetric building blocks. Here, we used 
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the symmetric XY-4 sequence as the building-blocks, and we modified the concatenation 
scheme in such a way that the symmetry is preserved and the delays between the 
pulses are identical at all levels of concatenation. The conventional (asymmetric) version 
CDD„(A) is iterated as [CDDn-i-X -CDDn^i-Y]"^. In contrast to that, we construct 
the symmetrized version CDD„(S) as yCDDn-i - X - CDDn-i -Y - y/CDDn-i\ ^ [22]. 
For n = 1, we have CDDi = XY-4. In Fig. [10] we compare the process fidelities |l5l |16] 
for the two versions of the CDD-2 sequence. Clearly, the symmetrized version CDD2(S) 
shows a significantly improved performance, compared to the standard, asymmetric 
version CDD2(A). 
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Figure 10. Process fidelities of tlie symmetric vs. asymmetric version of CDD-2 as a 
function of time for an average delay t — 12.5/is between the refocusing pulses. 

The results presented in this paper show clearly that it is time reversal symmetry 
is a useful tool for improving dynamic decoupling sequences. The symmetric sequences 
often perform better and never worse than non-symmetric sequences, at no additional 
cost. 
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Appendix A. Average Hamiltonian Calculation 

In this Appendix we sliow liow tlie average Hamiltonian was calculated. The essence 
of average-Hamiltonian theory is that a cyclic evolution, U{t), can be described by 
an effective evolution governed by a time-independent Hamiltonian H. When U is 
comprised by a sequence of unitary operations, i.e. U = e^^ ■ ■ ■ e^'^e^^ , the average 
Hamiltonian can be approximately computed by recursive applications of the Baker- 
Campbell-Hausdorff formula 



To approximate the average Hamiltonian of the XY-4 sequence, consider the 




i([A, [A,B]] + [[A,B],B]) 
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following sequence: 

[Ti-Ri-r-R2-T-R3-T-R4-Tf], (A.2) 

where Ri = R3 = R{X), R2 = R4 = R(Y) and -R(0) is a the pulse propagator defined 
in Eq. [3l In the asymmetric form of XY-4, Tj = r and Tf = 0. For Tj = rj = r/2, we 
obtain the symmetric form of XY-4. The total sequence propagator is 

U = e-'^^f (H Rke"'^A Rie~'"^\ (A.3) 

\fc=2 / 

where H is the Hamiltonian of Eq. ([T]) . 

To account for fiip angle errors, we decompose the propagator into a product of the 
ideal pulse propagator sandwiched between two additional evolutions: 

R^ = e-^(^+^)"^^ 

= Q-iH/-^^-i7vS4,^-iH/-§- (A. 4) 



where = f^S^ and tp is pulse length. Substituting equation flA.4p in flA.SP and using 
the following approximation: 

^-iH/-^^-iHT _ ^-i{HT+H^*-§-) 

~ g^-i{H+^e'rrS^)T 



\ (A.5) 
the new sequence propagator is then rewritten as 

U = e-'^'^^f (n RkC-'^'A Rie~'^'^^\ (A.6) 



\k=2 



where 



nH[ = T,H + (A.7) 

^^^1=2,3,4 = rH+^jiHx + Hy) (A.8) 

TfH', = nH + ^-^Hy. (A.9) 

The calculation can be simplified by transforming the Hamiltonians to a new frame after 

each pulse, the so-called toggling frame. The Hamiltonians in this new frame are 
given by: 

rOh =n{HE + HsE) + \Hx (A.IO) 
rHk=2,:,A = r[HE + {-1)''+'Hse] 

- tp[(24,2 + l)Hx + (25fc,4 + l)Hy] (A.ll) 

TfH, = TfiHE + Hse) + |i^y (A.12) 
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and the sequence propagator is: 



5 



f/ ^ n 6"^^'=^'= (A. 13) 

k=l 

The final step consists of the recursive apphcations of eq. f lA.l|) to eq. (]A.13|) . An 
exphcit calculation of the zeroth and first order terms leads to the equations (jlj) and 
([5]) for the symmetric case and (l6l) and (I7l) for the asymmetric case. The calculation for 
XY-8 follows the same procedure as described for XY-4. Here the sequence is comprised 
by eight pulses and nine delays (see table [1]), this leads to the total propagator analog 
tolAlel 

U = e-'^'^^f ({[ Rke~'^'A Rie~'"'^^\ (A. 14) 



Vfc=2 



Transforming the Hamiltonians H[. to the toggling frame, one can show that Hk = Hi^k, 
the Hamiltonians i?fc=i,2,3,4 are the same as in Eqs. (lA.lOp - Eqs. (]A.11|) and 
rlh = t{He + Hse) + tpHy. 



